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Abstract 

The thermodynamics of solvable isotropic chains with arbitrary spins is addressed by the 
recently developed quantum transfer matrix (QTM) approach. The set of nonlinear equations 
which exactly characterize the free energy is derived by respecting the physical excitations 
at T = 0, spinons and RSOS kinks. We argue the implication of the present formulation to 
spinon character formula of level k — 2S SU(2) WZWN model . 
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1 Introduction 



The ID spin systems have been providing problems of both physical and mathematical interests. 
Among them, there exists a family of solvable models of Heisenberg type with spin-S" ]l|, . This 
includes, for instance, 
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as S = 1/2, 1, respectively. 

Ground state properties as well as low lying excitations have been elucidated by the powerful 
machinery of solvable models, the Bethe Ansatz Equation (BAE). It has been demonstrated in 
many contexts §§ |, | that the underlying field theory is level k = 2S SU(2) WZWN model §. 

Although this ID quantum model is equivalent to a 2D vertex model, physical excitations 
carry both natures of vertex models and restricted SOS models. This is firstly demonstrated in 

based on the S— matrix argument. The space of states is identified in Q. By decomposition 
of crystals, such double feature is made explicit in terms of "type of domain walls" and "type 
of domain". There are independent justifications of this: The double feature in the spectral 
decomposition is shown by path space approach Q . See also fl(| for the decomposition of space 
picture realized in Fermionic forms. 

Here we are interested in the finite temperature problem. Standard arguments employ the 
string hypothesis^, [Tlf| . The excitation is described not in terms of "physical excitations" in the 
above sense, but in the "string basis" . There strings of arbitrary lengths are allowed, which results 
into infinitely many coupled integral equations among infinitely many unknown functions. The 
description successfully reproduces the expected specific heat anomaly. It may not be, however, 
best suitable for practical numerics. 

We revisit the problem via the recently developed commuting Quantum Transfer Matrix 



(QTM) approach [13|- |20f| . The formulation does not reply on the string hypothesis. Rather, 



it only relies on analyticity structures of the object called QTM[21, IT^]. The problem of the 



combinatorial summation, i.e., evaluation of the partition function, then reduces to investigations 
on the analyticity of suitably chosen auxiliary functions. Up to now, two kinds of choices are 
adopted independently. 

(A) The eigenvalue of the QTM as given by the Quantum Inverse Scattering Method consists 
of several terms. The auxiliary functions are chosen from combinations of products of these 



terms |13| , 14, 15, [17]] . A convenient choice leads to a finite number of coupled non- linear integral 
equations for a finite number of unknown functions 

(B) A set of auxiliary functions may be chosen from the fusion hierarchy among "generalized" 



QTMs [12, 16, Il8|, [2(1. Generically, one obtains an infinite number of coupled non- linear inte- 



gral equations for an infinite number of unknown functions. This can be shown to recover the 
conventional TBAs based on the string hypothesis. Of course, the new approach is entirely free 
of any assumption about excitations like string hypothesis. 

The spin-1 case is analyzed in the related problem, in the context of finite size corrections p3| . 
There six functions are introduced in the spirit of (A). The structure of NLIE among them is 
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much more involved in comparison to the spin-1/2 case, and seems to defy a simple-minded 
generalization to higher S cases. In a sense the most subtle point in the QTM approach appears; 
one does not know a priori "better" set of auxiliary functions. 

Here, we adopt other choice of auxiliary functions, in particular for the spin 1 case the 
number of these functions is 3 in contrast to 6 as in [p3| . A simple idea of combining the two 
formulations (A) and (B) works well so that the generalization to arbitrary S is possible. The 
adopted functions agree with the picture in Q. Roughly speaking, the fusion part (B) of the 
auxiliary functions is related to the RSOS piece of the excitation, while auxiliary functions from 
(A) correspond to spinons. 

We remark that the fusion hierarchy itself is not truncated, by brute force, into finite set. 
Instead, the spinon part makes the functional relations among them strictly closed. Thus we 
obtain 25* + 1 coupled integral equations for 2S + 1 unknown functions. 

Besides the practical advantage, it implies the universality in the description of thermo- 
dynamics of solvable quantum ID chains, i.e., the description only in terms of objects which 
reduce to physical excitations in T — > 0. This has been already demonstrated for several models 
in highly correlated ID electron systems including the supersymmetric t — J model [14|, the 
supersymmetric extended Hubbard model [^] and the Hubbard model [17]. There the exact 
thermodynamics are formulated in terms of "spinons" and "holons", although they lose sense 
at sufficiently high temperatures. The present study adds one successful example even in the 
fusion models and gives further supports on the above conjecture. 

This paper is organized as follows. In the next section, we define the main object in this 
approach, the quantum transfer matrix (QTM). A minimal information of the novel approach is 
sketched. Section 3 is devoted to a brief description of the fusion hierarchy of generalized QTMs. 
After these preparations, we introduce auxiliary functions and examine functional relations 
among them in Section 4. The analytic structure studied numerically leads to nonlinear integral 
equations as discussed in Section 5. Based on these equations the low temperature asymptotics is 
studied analytically in Section 6. The central charge of the level k = 2S SU(2) WZWN model is 
successfully recovered. We also present the numerical evaluation of specific heat of S = 1/2, 1, 3/2 
models for wider ranges of temperatures. In Section 7, the implication of the present formulation 
to the spinon character formula of the WZWN model [24, ||, |4^[ is discussed. The summary of 
the paper is given in Section 8. 



2 QTM Formulation 

The present QTM formulation originates from two ingredients: the equivalence theorem between 
ID quantum and 2D classical systems [^] on one side and the integrability structure on the 
other ^5|. Especially, the latter provides the way to introduce commuting QTMs which reduce 
the problem of combinatorial counting to that of analyticity of suitable auxiliary functions^]. 



Such a strategy has been successfully applied to several interesting models ]lj]-[^0|. We also 
mention earlier studies on thermodynamics |26|] -|32| which essentially utilize only the former 
part of ideas. 

A classical analogue to solvable spin-S" XXX model is already found as a 2S + 1 state vertex 
model The Boltzmann weights are identified with the matrix elements of 5(2 invariant Ff 
matrix: 
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Figure 1: Graphical representation of the R matrix ( its element H^p 
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where k = 2S and i^- is the projector to Vj, the j + 1 dimensional irreducible module of sfa- 
We choose {—j/2, —j/2 + 1, ■ • ■ , j/2} as basis for Vj. The spectral parameter u represents the 
anisotropy of the vertex weights. The Yang-Baxter equation implies the commutation of row-to 
row transfer matrices for arbitrary spectral parameters u, v: T(u)T(v) = T(v)T{u) with 



^)=zn^3 + >)> w 

/i i=l 

where L denotes the real system size, aj, /3j, m G 7Z = PR y and P(x ® y) = y <8> x. 
The Hamiltonian is obtained as the logarithmic derivative at it = 0, 

n = j— inr(«) . (5) 

du u=o 

It is an easy exercise to verify @ gives (Q) and( |2|) for S = 1/2,1 respectively. This may be 
most easily done by representing Pj in (||) by 

ft nn 

p=0,p^j 

and = 1/2(|(|-|-1) — A;(| + l)). The Hamiltonian for general S can be extracted similarly. This 
is the well-known expression of the equivalence between ID quantum systems and 2D classical 
models. To utilize the equivalence in evaluating finite T quantities, especially, free energy, we 
need to proceed further. 

Let us introduce Boltzmann weights 1Z (1Z ) of models related to (||) by clockwise (anticlock- 
wise) 90° rotations 
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Figure 2: In the same spirit to FigjT], R matrices, 1Z and 1Z are depicted as above. 



The standard initial condition of the i? v matrix and (||) lead to significant relations, 



T(u) = T R e 



uH/J+0(u 2 ) 



T(u)=T L e 



uH/J+0{u 2 ) 



(6) 



where T is defined in analogy to (|J) , replacing 1Z by 1Z. Tr,l are the right- and left-shift 
operators, respectively and they commute with the Hamiltonian. 
We are ready to apply Trotter formula; by substitution 



u = -J/3/N, 



we find 



< — \N/2 

T{u)T(uY 



-m+o(i/N) 



(7) 



(8) 



where (3 denotes the inverse temperature. N is a large integer "Trotter" number, interpreted as 
a fictitious system size in a virtual direction. Thus the partition function of the quantum system 
(size L, inverse temperature (3) 



Z= lim Tr( T(u)T(u 



N/2 



(9) 



is identical to the partition function of an inhomogeneous 25"+ 1 vertex model with alternating 
rows on a virtual 2D lattice of size Lx N .See (Fig §). Although the above mapping is exact, 
the expression (||) is not yet efficient. The eigenvalues of T(u) T{u) are almost degenerate. Hence 
it is still difficult task to evaluate the trace. The intriguing point in is to consider a transfer 
matrix T(u) propagating in the "horizontal" direction. 

This novel operator acting on N sites, has gaps in the eigenvalues provided T > 0. Here we 



adopt further sophisticate approach developed in 12|- PQ| . Explicitly, we define QTM by 
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which reduces to above-mentioned T(u) by putting x = 0. Fig .|| represents QTM graphically . 
Though the introduction of the extra parameter x seems to be redundant, there is a remarkable 
property; 



[Tqtm{u,x),Tq TM {u,x')] = 
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Figure 3: Graphic representation of the partition function on a virtual two dimensional lattice 
of N x L. The operators T, T transfer states from bottom to top, while T(u) and Tqtm do from 
right to left. 
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Figure 4: Graphic representation of 7qtm(^>^)- 
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Figure 5: For Tj(u,x), vertical arrows carry Vfc while horizontal (broken) arrow does Vj. 



by fixing u. This originates from the fact that 7Z and 1Z operators possess the same intertwiner. 
Thus for each T, one can associate a auxiliary complex plane x to the partition function. 

Due to the gap in spectra, the free energy / of ID quantum spin chains is given only by the 
largest eigenvalue Aqtm(w, x), 
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This is the starting point of our analysis. The difficulty in evaluating (12) lies in the A" 
dependence of the vertex weights. The numerical extrapolation through finite A" studies may be 



plagued by marginal perturbations [19]. The prescription is to utilize the existence of complex 
plane x for each T. The analytic properties of "Tqtm and suitably chosen auxiliary functions in 
the x— plane make the evaluation possible and transparent. 

Before closing this section, we shall describe how to modify the above relations in the pres- 
ence of an external magnetic field H, namely by inclusion of the Zeeman term —2H Sf to 
Hamiltonian (||). This contribution is described by diagonal operator D(H), 



/ & -2Sf3H 
p -(2S-l)(3H 
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-2S/3H 
e -(2S-l)/3H 



V 
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Thus one only has to insert this inside the trace 



/ _ \ N/2 

Z= lim Tr[T(u)T(u)) D(H) 



In the rotated frame, the effect of the insertion of D(H) is translated to the boundary weight, 



N/2 

/i i=l 



U — IX). 



3 Fusion Hierarchy 

We consider a hierarchy of quantum transfer matrices acting on V® N . Let Tj(u, x) be a member 
of the hierarchy with the auxiliary space Vj . 
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In other words, it is the transfer matrix of the vertex model of which spins S(= k/2) are 
assigned to vertical edges and spins j/2 to horizontal edges (Fig||) The quantity of our interest, 
Tqtm(u,x) coincides with T)~(u, x) apart from over-all normalization, which is specified later. 

For brevity, we shall only give matrix elements of 7Z(u) defining the most fundamental 
T x (u,x). 

1lf) /2 ' ±l/2 (u) = u + l±2£, n^r/' e ~ e ' (u) = v / {k + 2 + 2mm{£,£'))(k-2min(£,£')) 



where £,£' € {—k/2 • • • , k/2} and \£ — £'\ = 1. Similar to (H), the corresponding R v matrix 
has decomposition, 

R v (u) = (u + k + l)P k+1 + {u-k- (13) 

TZ— matrix for Tj(u, x) is obtained from the above elementary TZ(u) by j — 1 times fusion in the 
auxiliary space. By the construction, arbitrary pairs in this hierarchy are commutative if they 
share the same u, 

[Tj(u, x), Tji(u, x')\ = 0. 

This is a generalization of (|ll]) . In the following, we fix u for all QTMs and omit the dependency 
on u. Due to the consequential commutativity, one needs not to distinguish operators Tj from 
their eigenvalues. 

Then the explicit eigenvalue of the most elementary transfer matrix T\ (x) reads 
Tx(x) = - (k - - (k + fy)e? H Q ^Q(J l) 

+4>-(x + (k - l)i)Mx + (k + l)i)e-^ H Q{ - X ~^ 

N/2 



(j)±(x) := (x±iu) 
m 

Q(x) := JJ(x-Xj) 

3=1 



(14) 



where Xj, (j = 1, • ■ • , m) denotes the solution to the Bethe ansatz equation: 

(f>-(xj + (k - l)i)<f>+(xj + (k + l)i) _ 2i3HQ( x j + 2 
4>+(xj — (k — l)i)(j)-(xj — (k + Q(xj — 2i) 

The number of BAE roots, m, differs for different eigenstates generally, and m = Nk/2 for the 
largest eigenvalue case. 

By construction of the fusion hierarchy and from singularities of the intertwining operators 
(|l~3| ) by tentatively replacing k — > j, the following relation is valid: 



Tj(x)Ti(x + 1)) = Tj +1 (x-i)+gj(x)Tj-i(x + i) 

gj (x) = 4>_(x- i(k + j + 2))4>„(x + i(k - j))4>+{x - i{j + k))(f> + (x + i(k -j + 2)). 



S 



Prom this, one can prove the important functional relation (T— system) by induction |3 

T p (x + i)T p (x - i) = f p (x)+T p - 1 (x)T p+1 (x), (p>l) 

p 

f P (x) ■= J] II Mx + Mp- k-2j + l))(P a {x + ia(k- P + 2j + 1)) 

3=1 c=± 

(15) 

and To = 1. 

By substituting ( |14] ) into ([H]), we can successively obtain T p (x),(p > 2). Explicitly, T p (x) 
consists of a sum of p + 1 terms, 

p+i 

T p (x) := J>f(x), 

£=1 

y,<P) M . = r WiP+2-2l) ,(p) M Qix + i(p+l))Q(x - i(p+l)) 
e W • Wi [ iQix + i^-p-l^Qix + i^t-p-?,))' 

p-e+i 

i>t\x) := Yl (t>-{x + i{p-k-2j))(j )+ {x + i{p-k + 2-2j)) 

i=i 
e-i 

xY\4>-(x - i(p - k + 2 - 2j))<p+(x -i(p-k- 2j)). 

(16) 

As is previously noted, Tk{x) has a normalization trivially different from Tqtm(^ ; ^) in the 
previous section: 

Tfc(x) 

rQ ™ (n ' x) = Y^M 

fa(x) ■■= * N/2 

(17) 

In the original problem of the spin-S 1 chain, only T^ix) is of our interest. The auxiliary Ij's, 
however, make the evaluation of Tk{x) transparent as is shown in the following. 

4 Auxiliary Functions and Functional Relations among them 

To explore the analyticity of the transfer matrix T^(x), we introduce k + 1 auxiliary functions. 
The first k — 1, functions, {yj(x)} are already found in several literatures and have sound basis 
in the sfa fusion hierarchy. They are defined by [36, [37]] 



wW = jrw ■ > £ 

The functional relations among them are sometimes referred to as Y — system: 

yj(x + i)yj(x-i) = Yj-i(x)Y j+1 (x), j>l 

Yj(x) := l + yj (x) (18) 
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and yo(x) := which is a consequence of eg. (|15D . Note that the Y — system is not truncated to 
a finite set in this case. The k — 1-th equation, which characterizes yk-i{x) inevitably contains 
Uk(x) in r.h.s., and so on. Thus another device is needed to construct a finite set of auxiliary 
functions satisfying a complete and closed set of functional relations. The remaining two func- 
tions b(x),b(x) and their "relatives" 23 (x) := 1 + b(x), 23 (x) := l + b(x) play this role. We define 
them by ratios of A's in T)%(x) as, 

v _ xP(x + i) + -.. + \P(x + i) 
[X) '~ \( fc ) i , ^ 



b(x) 



xi k \x- i ) + ... + xi%( 



\[ k \x-i) 



(19) 



The following relations are direct consequences of the above definitions: 

BOOA&C* + = e-W H <B(x) [J fl + (2J + *)0 %£~)5 = T ^ x + 

<7=±j = l ^ ' 

iJ L AJ : Q(x-tk) 
cr=± j=l v 7 



(20) 



We have fe — 1 equations for y^, (j = 1, • • • , k — 1) in terms of Yj(x), (j = 1, ■ ■ ■ , k — 1), 23(x) 
and 23 (x). The first k — 2 equations are chosen directly from the Y— system. In the k — 1 th 
equation, we rewrite Ifc(^) in the rhs of Y — system (j = & — 1 in ( [Tq) ) by 23(x)23(x), thanks to 
( pOj ) , the definitions of y^ , Y^ and the functional relation ([D]) : 



y fc _i(x - i)y*-i(a; + i) = r fe _ 2 (x)23(x)23(x). (21) 

Finally, equations for b's in terms of Yj(x), (j = 1, • • • , k— 1), 23 (x) and 23 (x) are to be found. 
By comparing explicit forms, one finds 

b(x) _ 0(k+i)H TT Mx + m) Q(x + i(fc + 2)) 

b(x) = e -WTT Q(x-^ + 2)) 

Il* =1 CT (x - (2j - a)i) Q(x + l k) 

(22) 

Note that r^_i(x) is presented by Yfe_i(x): 

T fe _i (x - i)T fc _i (x + i ) = /fc-i (x) Yfc_i (x) (23) 
which originates directly from definitions of yk-i,Yk—i and the functional relation fll5|). 
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In what follows, we analyze these functional relations via the Fourier transformation. One 
denotes d/b[g] to mean the Fourier transformation of the logarithmic derivative of b(x), 

dlogb(x) Jqx 



dlb[q] := / ^^e lqx dx, 

and similarly for other functions. Under some assumptions on analytic properties of the auxiliary 
functions, the above functional relations can be transformed into algebraic equations in the 
Fourier space. Roughly speaking, one can solve dlQ[q] functions in terms of dl^&[q] and dl*B[q] 
by deleting dlTk[q] from algebraic equations originated from (|20|) . Similarly dlTk-i[q] is solved 
by dlYj e _i[q] from (|23|). Substituting these results into Fourier transformations of logarithmic 
derivatives of (f22|), one finds expressions of dlb[q], dlb[q] in terms of dl*B[q], dl^[q\ and dlYk-\[q]. 
After the inverse Fourier transforming and the integration over x, we shall find the desired finite 
set of equations. We will make the above mentioned analytic assumptions explicit and examine 
them in the next section. 

Before going into details, let us discuss the physical interpretation of the above functions. 
As is argued in 0, the S— matrix of excitations in the spin-5 model factorizes into two pieces: 
spin-^ SU(2) S— matrix and RSOS S— matrix of s/2 level k = 2S . This is consistent with the 
general expectation that the underlying field theory is the level-A; SU(2) WZWN model. The 



latter is known to "decompose" into Gaussian and Z^ parafermionic field theories [34]. One 



finds, see for instance [35], several evidences for the equivalence between the SI2 RSOS model 
in regime II and the Z^ parafermion field theory in the scaling limit. In the present description, 
b, b are to be identified with up- and down- spinons. As we will see later, only they couple to 
the magnetic field directly. For the S = 1/2 case, there are further direct evidences for this 
identification |5l]] . On the other hand, {yj{x)} are insensitive to the external field. We are led 
to identify yj{x) in our choice as the RSOS piece of the excitations. The RSOS model possesses 
a subset of the Y— system fll8|). The additional condition = for the model leads to the 



truncated set of equations among k — 1 y's |38|, pq| . In the present problem, at sufficiently low 



temperatures, one observes, |b|,|b| ~ for x « ln/3. Thus the substitution of 23 = 23 = 1 
in (|2l| ) might be legitimate in the vicinity of the origin. Then the resultant approximated Y — 



system coincides with that of the RSOS model. In this sense, (21) represents a gluing relation 
between spinon and RSOS parts of excitations. 



5 Nonlinear Integral Equations 

We derive the nonlinear integral equations among auxiliary functions introduced in the previous 
section. The crucial observation is, that all nontrivial zeros and singularities of these functions 
are determined by zeros of Q(x) and Tj(x),(j = 1, ■ • ■ ,k). For the largest eigenvalue sector 
of Ifc(x), zeros of Q(x) form so-called k— strings. Imaginary parts of zeros are approximately 
located at (k + 1) — 2i, £ = 1, • • • , k. For later use, we introduce notations, 

:= Q(x-ik), V 2 (x) ■= Q{x + ik). (24) 

Empirically, similar patterns are found for zeros of Tj(x): they distribute approximately on lines, 
Qx = ±(k + j — 2£),£ = 0, • • • , j — 1. We assume that these observation from numerics with 
fixed N is valid and that the deviations from lines are very small in the limit — > 00. Then one 



11 



deduces the following ansatz on the strips where auxiliary functions are Analytic, NonZero and 
have Constant asymptotic behavior (ANZC). 



b(x),«8(x) - 1< 3z < 0, 
b(x),tb(x) 0<3x<l, 

y j (x),Y j (x)(j = !,■■■ ,k-l), T p (x), (p = 1, ■ • • , k) - 1 < Gx < 1 
*i(x), 3=x < 
i ff 2 (x), 3=x > 

We find it convenient to shift the definition of the arguments in b, 2$, b and 25. To avoid confusion, 
these new functions are denoted as o etc, 

o(x) := b(x — £7), o(x) := b(x + £7), 

and similarly for capital functions. Here < 7 < 1/2 is an arbitrary but fixed parameter. Note 
that this is equivalent to adopt small shifts in the definition of the integration contours for the 
Fourier transformation. Due to the ANZC properties of b, b, such modifications are almost trivial 
in the Fourier space. 

Having identified ANZC strips, we revisit eqs(pO|). Consider the integral, 



J^logT k (z)e igz dz, 

where C encircles the edges of "square" : [—00 — i, 00 — i] U [00 — i, 00 + i] U [00 + i, —00 + i]\J 
[—00 + i, —00 — i] in counterclockwise manner. Due to the ANZC property of 4- \ogT k (z) inside 
C, the following equation is valid from Cauchy's theorem: 

/OO J ,>oo J 

— logT fc (x - i)e iq( - x -^dx - — \ogT k (x + i)e iq( - x+i Ux. 
-00 ^x J— 00 

One substitutes eqs(pO|), rewritten in terms of ^^(x), 2l(x), 2t(x), into the above equation and 
derives identities among dl^ i^[q\, dl%[q], dl%[q\^ 



dl^ 1 [q<0] = 

(1- 
co 

,e~ kq sinh kq cosh(l + u)q 



^ P (l-7)9 _ P -(1-T)9 

dm[q>0) = —m[q]-- — — dm\ q ] 

I cosh q 2 cosh q 



+ iriN- 



cosh q sinh q 



e kq sinh kq coshfl + u)q 

— iriN 

cosh q sinh q 

m 2 [q>o] = 0, 
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and dl^i[q = 0] = —dl^>2[Q = 0] = irNki. Similarly, one can derive an identity for dlyj[q\s 
and dlYj[qYs from fllS|), and dlT k _i[q] and dlYf._i[q] from (^). Substituting these relations 
into the original definitions of a and o, we obtain k + 1 algebraic relations in Fourier space. 
(Remember ^1,2 (#) are related to Q(x) by After taking the inverse Fourier transformation 
and integrating over x once, we arrive at the k + 1 coupled nonlinear integral equations: 

/ Iogyi(x) \ / \ / logYi(x) \ 



log Vk-i{x) 
log a(x) 
V log o(x) / 





f3H + cf(u, x — i'y) 
\ -i9fr + d(u,s + i7) / 



logYfc_i(a;) 
logSl(aj) 
V log2t(x) J 



(25) 



where (/C * g)i denotes the matrix convolution, ^ ■ J* fCi t j(x — y)(g(y))jdy and the "driving" 
function d(u, x) reads 

A /" sinh uq 



d(u, x) 



-iqx 



dq. 



2 J q cosh q 

The integration constants (±/3H) are fixed by comparing asymptotic values (|x| 
sides. 

Explicitly the kernel matrix is given by 



oo) of both 



where 



/ 0, K(x), 0, 
K{x), 0, K(x), 

0, K(x), 

0, 
0, 
0, 

V o, 



0. 
0. 



o, 

K(x), 

o, 

0. 



K{x) :-- 
F(x) :-- 



0. 
0, 

0. 

K(x), 
0, 

K{x — i'y) 
K[x + i'y) 



1 



0. 



4 cosh ttx/2 
1 f 00 e~\ q \~ iqx 
2ir J-oo 2 cosh q 



0, 

K(x + i'y), 
Fix), 
-F(x-2i(i- 7 )), 





if (a; — £7) 
+ 2t(l - 7)) 
F(x) 
(26) 



Tfc, in terms of these auxiliary functions, can be derived similarly. Technically, we find it conve- 
nient to introduce, 

T k (x) 



T k R (x) :- 



(27) 



Y\ p=1 (t)e^){x - 2i(k + 1 -p))<j> e (p + i){x + 2i{k + l- p)) 

and e(p) = +(— ) for p = even (odd). Then the product of two equations in ( |20| ) leads to a simple 
algebraic relation: 



T?(x-i)T?(x + i 



j *B(x)*B(x) for k even 



»oo»(*) £K%:fcij for odd 



(28) 
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Figure 6: specific heats for S = 3/2, 1 and 1/2 from top to bottom 



From ANZC property of the both sides of (28) in appropriate strips, the logarithmic derivative of 
them reduces to a simple algebraic equation in Fourier space. Taking account of the normalization 
(H7D and (|27h , we have 



\ogk QTM (u,x) = \ogh^ TM {u,x) 

+ r loggt(y) + r loggtfo) 

7 4cosh7r/2(x — y + 27) y 4cosh7r/2(x — y — £7) 

logA^ M (u,x) := -J fc=1(mod2) ^ J e -^-M^±d q (29) 

* (/> e(p) (x - 2i(k + 1 - p))<ft e (p+i) (x + 2i(fc + 1 - p)) 

P 1 

(30) 

Finally put u = —(5J/N and send A" — » 00 analytically. This merely amounts to replacements: 

/?Jvr 



d(u,x) — > D(x) 



2cosli7r/2x 



logA§r M (u,x = 0) -» -/3e = — — ^ ; + ^ = i {mod2 )/3«/log 2. (31) 

i=i 3 

Note that eo coincides with the known ground state energy after a trivial shift which stems 
from the difference in the normalization of R v . The k + 1 coupled nonlinear integral equations 
and log Aqtm do not carry the fictitious parameter N any longer. They efficiently describe the 
thermodynamics of the solvable spin-S" XXX model. For an illustration, the specific heat for 
S = 1/2,1,3/2 is evaluated for a wide range of temperature and plotted below. Each curve is 
produced by 10-30 minutes CPU time calculation on a Micro Sparc Work Station. In the next 
section, we derive the low temperature properties using (25), (|29| ) and (|3ll). 
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6 Analytic Evaluation of the Low Temperature Asymptotics 

We consider T — > for the vanishing magnetic field. In sufficiently low temperature regime, 
a, log 21 shows a sharp crossover behavior like a step function: |a| , | log 2t| <C 1 for \x\ < ^ log 7T/3J 
and |a|, | log 2t| = 0(1) for \x\ > - log7r/3J. Thus the following scaling functions [23] control the 
asymptotic behavior: 



loga(±-(£ + log7r/?J)), 

7T 

iogo(±-(e + io g ^j)), 

7T 

2 



L4±(0 :=log2l(±-(e + log7r/3J)), 

7T 

ZA±(£) :=log2t(±-(C + log7r/3J)), 

7T 

2 



Zy±(£) := logy p (±-(£ + log7r/3J)), IY±(£) := lo g y p (±-(£ + logic /3 J)). (32) 

7T 7T 



In terms of these scaling functions, NLIE are expressed by, 



I 











_ e -£±*77i-/ 2 
\ _ e -?T«7T/2 

IT 7T 



ia±(o 
V ^(0 J 



(33) 



Note that neglect of small corrections ~ 0(T) leads to the decoupling equations for ±. 

The thermal contribution (the second and the third term in (p9|)) to limjv^oo log Aqtm(^ 



■jj, x) reads, 



„nx/2 r [■ 

l__[ e i^/2 j e-ZlA+dt + e-^ 12 J e'^lA+d^] 
+^7[ e " i77r/2 y e^lA-di + e 1 ^' 2 J e-HA-di]. 



(34) 



The crucial observation in [23] is that one needs not solve (|33| ) to evaluate (34) provided that 
the kernel matrix function satisfies a symmetry, /Qj(x — y)= fcj,i{y — x). This property is valid 
in the present case. See (|T~ 
We define F± by 



00 k 1 d 



d 



-°°p=i ^ ^ 



+ / [(^-la^lA* + Aa^lA^ - (i^)^ - ^-lA^la^. (35) 
'-00 <K <K «? 



Then the trick in [|23|] is as follows. First, take the derivative of both sides of (|33|) and multiply 
them by a row vector 



(36) 
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We call the resultant equality (A). Second, multiply both sides of ( |33| ) by the derivative of the 
row vector (|36|), which is referred to as (B). Finally, subtract both sides of (A) and (B), and 
integrate over £. Then the lhs of the resultant equality is nothing but F±. Remarkably, the most 
complicated terms in the rhs, like 

and 

cancel with each other. To be precise, the first integral can be converted step by step, 

= -JdtdeiY+(o d ^~ i Y ;(o 
= f dtdt'iY+m^i '-t)±tY+tf) 
= J d^^iY+m^-ew+io 



where the symmetry of the kernel matrix, partial integration and the change of integration 
variables £ <-> £' are used. 

Similar cancellation happens for other terms, and the following equality results, 

F± = 2 J \e~t ±i ^l' 1 lA ± + e-teWHA^dt. (37) 

The first thermal correction is thus given by 

e irx/2 p + e -nx/2 p 

~2^T + 2^(3 J • (38) 

To evaluate F± explicitly, one rewrites the integration variable from £ to a,a,y p . For example, 
the first summation term in F± is transformed to 

r p±(00) a A°s(l + y) logy^ v\r r , ±, ^ T ( ±( m 

2^ / + d v\ — : 7—7) = 2 z^{ L +\y v (°°)) - L +(y P (-°°))}> 

p =i ^ (-00) y L ^ y p =i 

and similarly for others. L + (x) is a dilogarithm function and is related to Rogers' dilogarithm 
function L(x) by L + (x) = L(x/(1 +x)), 



L+(x) :-- 
L(x) :-- 



1 f x log(l + y) logy,, 

V) + ^)dy 



2./0 ' y 



1 

'2. /n 



-logfl 



(39) 
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The asymptotic values of scaling functions are easily extracted. For x — ► oo, coincides with 
original b. Thus one derives the limiting value by its definition ( |i~9| ) in terms of Ap 's. Similar 
for y^(oo). For x - 



,„ . ... , -oo, one should rather consult (33|). We send the argument x ■ 

both side and solve the resultant algebraic equations. The results are summarized as, 



-oo m 



-oo) 
-oo) 



o ± (-oo) = 0, 

np ■ ir(p+2) 

sm ^ sin 



8111 



a ± (oo) = a ± (oo) 
1 <p < k-1, 



k. 



k+2 

p(p + 2) l<p<k-l. 



(40) 



With these pieces of information, F± is now given by 

k-1 



P =i 



Finally we use three relations [| 



Pip + 2) 
(p + 1) 2 



Li 



sin JIE. s i n ^(p+ 2 ) 
bnl k+2 blIi k+2 



sm 



2 7r(p+l) 
k+2 



)] + 4L( 



1 + 



L(l) 



3n 



n + 2 

which yield the neat result 



L(x) + L(l-x) = — 
b 



x e [o, l] 



n-1 



j=0 

"Zl sin 2 - 



j=0 



sm 



2 7r(j+l) 
n+2 



7T 2 fc 



n € Z. 



>o 



Thus we conclude the low temperature asymptotics of free energy 



7T 



-c(fc), c(fc) = 3k/(k + 2), 



(41) 



where u s = Jtt/2 coincides with the known spin velocity[40] and c(k) is nothing but the central 
charge of level-A: SU(2) WZWN model. This is the desired result from the WZWN description of 
massless quantum spin chains. We remark that the final part of the calculation is quite parallel 



to that in [11] utilizing the same dilogarithm function identities. There the spin-S 1 XXZ model 
is discussed via the standard string approach at "root of unity" where the number of strings is 
truncated finitely from the beginning. 



7 Spinon Characters 

The character formulae obviously depend on the base of the space. Recently, the quasi-particle 
representation has been attracted much attention in the context of long-range interacting model 
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f_ 



Figure 7: A contour C[/_, /+|{2}|{-1, 1}]. 



[24, 41, |4^[ , spectral decomposition of path space in lattice models M , or in the statistical inter- 
acting picture of Bethe Ansatz solvable models p3|] -[47]. See also p8] , |49[ | in the different view 
points. For the spin-1/2 case, it has been shown recently that the novel thermodynamics formu- 
lation yields a natural spinon character and such character formula is generalized to sl(n)k=i 
[51]. Thus it is tempting to find analogues for sl(n)k=2S- The results given in the previous sec- 
tions provide the first step for the simplest n = 2 case as discussed below. Remark that we 
consider " chiral-half ' , such that the only + part contributions (say F + ) in the previous section 
are taken into account. 

The character needs the description of all excited states. In the present contexts, this infor- 
mation might be encoded in the additional zeros of a, a, y and their capitals in their "physical 
strips". Indeed, some low excitations are identified in such a way |i~6| , |Tg| ], and corresponding 
excited state TBA are derived. Such an analysis is of considerable interest, however it requires 
extensive numerical efforts. We leave it as a future problem and make a short-cut detour here 
employing the strategy in fl50| . 

The central charge is described by the dilogarithm function of which the integration contour 
is simple. On the other hand, one can define an analytically continued dilogarithm function 
Lq{z). This is established by adopting general contour C for the integration contour of the 
dilogarithm function. We then generalize a successful observation from specific examples; all 
excitation spectra, or effective central charge c e fr, in the conformal limit shall be described 
by Lc{z). Namely, the replacement of a simple contour in the integral representation of the 
dilogarithm function by complex one leads to an excited state. Regarding c e fj as a function of 
Lc(z), the summation of q~ c rff/ 24 over a certain set of C is thus expected to reproduce affine 
characters. (Readers should not confuse this formal variable q, the standard notation in this 
field, with the Fourier variable used in previous sections.) 

Let us be more precise. By C we denote a contour starting from /_ and terminating at 
/+, such that it crosses firstly [1, oo) 0) times then crosses (— oo,0] £i times then again 

r/2 times w.r.t. [l,oo) and so on. The intersections are counted as +1(— 1) if the contour goes 
across the cut [l,oo) in the counterclockwise (clockwise) manner and (— oo,0] in the clock- 
wise (counterclockwise) manner. (Note the definition is slightly different from p0| , |5T[.) We 
denote this by C[/_, £2) • • • }|{?7lj ?72> • • • }]■ The set of contours are parameterized by 

S = {Cif^jM ,\V?,--- ,]}, where = y+(±oo)/Y+(±oo), (1 < p < k - 1), 

= a+(oo)/2l+(oo), /| fc+1) = a+(oo)/2t+(oo) and / W = f[ k+1) = 0. 
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In the absence of additional zeros of auxiliary functions in the " physical strip" , we have 

= log(/| p) ) - £ 0py log(l - fW), (42) 

where the "statistical interaction" matrix q is related to the zero mode of the Fourier transfor- 
mation of the kernel matrix (|26|) by 

fl = /-/C[g = 0]. (43) 



This is the situation we have treated in previous sections, and (42) follows from (p3|). One 
replaces log in (f42|) by an analytically continued one, Log c (z) in excited states: 

nD® = Log c(p) (/| p) ) - ^ 0py Log c(pO (l - /| p,) ), (44) 

p> 

where Log C ( P ) (/if ) = log(/+^) — 2niz*2t€i and so on. Here is introduced for consistency 
of both sides, and is interpreted as "chemical potential" On the other hand, should 
originate from the zeros of auxiliary functions in the physical strips. One may be able to prove 
that such actually agree with one in j44|), in principle. Though such microscopic derivation 
is yet to be done, we assume the coincidence in the following. 

Now, the excitation spectrum is solely implemented in the effective central charge c e g(S) 

c cff (5) = ^£(W/^/? } ) - y^ (p) Log c(p) (l - ff)). (45) 
p 

Here I)W term is included by hand, so as to match its interpretation as chemical potential [pQ], 
L c is given by, 

icU-M = -5 / {^c^A + ^M )dz . (46) 

z j c Z 1 — z 

for a path C from /_ to f + . 

After straightforward manipulations, one finds |5 



L CW (f M , f?) = L(f®) - L(/«) - ™ ^ d P) log(l - /?) - vri £ V? log(/f) 

+2- 2 (E d p) x£ ^ - 4^ + • • • + 4 a) y 

ii i 

(47) 



Remember L(z) is defined in (^). The substitution of (|47| ) in (flq), using ( f44| ) and the definitions 
of Log C ( P ) leads to a remarkable result: c c g- can only be written in terms of c(k) (see ([41])), } 
and {r^ 1 }, 



fe+1 

c cfr (5) = c(k)-24T(S), T(5) = -*n n + EE^ a) ( r ?i a) + --- + ^ a)A 

a=l £>1 

*n = (n (1) ,-- - ,n (fc+1) ) (48) 
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Figure 8: Young diagrams YD^ of = 3 corresponding to a set of winding numbers: {£} = 
{1, 1, 1}, M = {1, 1, 2} (left), {£} = {1, 1, 1}, {r/} = {0, 1, 1} (middle) and {£} = {1, 1, 1}, {r,} = 
{1,1,0} (right). Arrows are put for a guide to eyes. The most right diagram is termed "with 
tail" in the text. 



where = rffi . Note that the explicit forms of are not needed in the above trans- 
formation. 

In the following, we shall argue that the summation of q~ Ccii over some subset O of all 
possible contours reproduce the character ch j(z, q) of level k WZWN model with spin-j (j = 
some fixed integer or half- integer) . 

We present necessary conditions for such O below. 

• rj's and £'s are non-negative. 
Such path can be parameterized by 




and^ p) >0,(p = l,--- ,nW). 

• For p = k, k + 1, we require " (fc) +^ (fc 1} _j g z>o in addition. 

. For p < k - 1, C% > 1- 

Graphically, one can associate a Young diagram to a set of winding numbers: }, {Ve } 

(or a Young diagram with tail for some cases in p = k, k + 1). First draw a line of length 

(p) 

1 downwards. Next draw a line of length ^ to the left. Then again draw a line of length 
1 downwards. We continue this procedure times. Finally draw a horizontal line from the 
starting point to the left and also draw a vertical line from the end point upwards. See the Figure 
||. Obviously, the number of boxes in the diagram is equal to the second term in T(S). 

We allow for contours which are isomorphic to a set of Young diagrams YDb\(p= l,---k- 
1) such that 

s n (k-2p+l) _ n (k-2p) _ even _ 

• By the definition, the depth of YD^ is tv?' . The width is restricted by the maximum 

(p) 

value wiriax which is determined by depths of "adjacent" diagrams: 



,00 

7 max 



, v = l/2(n (p -^ + - 2n® + <%, fc _ p ), (p =!,-■■ ,k 



for the fixed j and = . The case p = k — 1 is exceptional, i 1 ' = l/2(n^ k 2 ) + 
n (k) + n (k+i) _ 2n (fc-i) + S jtl / 2 ). 
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Under the above restrictions on O, we find 



-A(i)+;/2 chj(Z;(?) 



^ 9 -c cff (5)/24 z (nW-n( fc + 1 ))/2 



seo 



n {h) n (k+i) >q 



SM,N(z;q) stands for contributions from (M,N) spinons; 

1 



S M ,N{z;q) 



{q)M{q)N 



(M-N)/2 



resulting from summations over "nodes" > 0,1 = 1, • • • and p 

Cartan matrix for slk+i and ^ A k denotes 



k, k + 1. Ai. is the 



l/4m.Aj.m 



q 

ra2,m3,- ,m fe i=2 



n 



|((2 - Ajfc).m + Uj)i 



mi 



and (itj)i = <5i,2j+i- The summations are taken over odd (even) positive integers for 



^even(odd) ■ 



Note we redefine nA fe ^ = mi+i, 1 < ^ < — 1. The appearance of the Gaus- 



sian q-binomial usually originates from combinatorics on the truncated Bratteli diagram and is 
the reminiscence of the RSOS model. Here the origin is also simple. It comes from the restriction 
on the width of Young diagrams. We denote the number of boxes in a Young diagram 
by 6(f). For fixed {nW}, T(S) assumes the same value for diagrams having identical b^\ This 
multiplicity is given by p(n^ p \w^l x ,b^) = the number of partition of into at most 
part, each < n^ p \ Thanks to the generating relation, 



(p) 



6b) 



If"' + Wmax 
n (p) 



In 



(p-i) 



n 



(p+i) 
„(p) 



we obtain the product of Gaussian q-binomial equivalent to the one in \&™* (m; g). 

Our result implies that the Hilbert space of level SU(2) WZWN model is isomorphic to a 
part of homotopy space of Rogers' dilogarithm function. We believe that this provides a novel 
interesting view point in the prototype of CFT and deserves further examinations. 



8 Summary and Discussion 

In this report, we formulate a novel description of the thermodynamics of solvable spin-S 1 XXX 
models. The suitable choices of auxiliary functions yield a natural generalization of the strat- 



egy in [ 23 1 . The nonlinear integral equations close finitely, which clearly differs from the string 
formulation 0] and is obviously efficient in numerics. The resultant formulation has an interpre- 
tation in terms of physical excitations, spinons and RSOS kinks. This has been demonstrated 
by the calculation of the low-temperature specific heat as well as the spinon character formula. 
The latter, however, is derived under several assumptions on the possible homotopy class of Lq- 
Certainly these constraints on winding numbers have "microscopic" origins in patterns of zeros 
of y— (or T) functions. This has been actually demonstrated for few cases: the superintegrable 
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3 state chiral Potts model |)2| and the simplest case of sl% RSOS models with open boundaries 
1 53]. We hope to report on extensive numerical investigations on zeros in the present context in 
the near future. 

Finally we mention spinon pictures arisen from different view points [54, 55, 56]. where 
explicit "spinon" bases are constructed by vertex operators. Thus the meaning of "spinon" is 
more transparent in comparison to the present approach. Less obvious is their concrete relation 
to eigenstates of spin Hamiltonians. The method also involves an uncontrolled approximation, 
a truncation procedure, in evaluating the partition function as well as "one body" distribution 
functions. This contrasts to the present approach which involves no approximation. In a sense 
they are complementary, and their relations are to be explored. 
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